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I . INTRODUCTION 

In  1977,  Lloyd  and  Brown1  investigated  the  feasibility  of 
controlling  a 105mm  spinning  projectile  by  means  of  horizontal  and 
vertical  side  forces.  Their  numerical  calculations  yielded  the  sur- 
prising result  that  an  applied  constant -amplitude  yaw  moment  could 
cause  dynamic  instability.  The  usual  linear  analysis  seems  to  predict 
that  such  a moment  would  cause  a steady-state  horizontal  trim  angle 
but  would  nave  no  effect  on  the  dynamic  stability. 

vatioIhtL,frfhfiC;j-^  wai\refolved  and  Brown  through  the  obser- 

ation  that  the  differential  equation  for  the  angular  motion  in  fixed- 

P ane  coordinates  contained  nonlinear  terms  in  Jpp  , the  coordinate 

!^-SPin  r^e'  These  usually  neglected  terms  vanish  completely 
when  the  equations  are  transformed  to  nonrolling  coordinates.  The 

non l*i near °terrnf  IJV horiz?,-tal  fnd  vertical  control  moments  become 
nonlinear  terms  that  can  be  easily  linearized.  The  resulting  sixth- 

flcZrTs  > aPPr^iniately  solved  for  large  gyroscopic  stability 

tainpd  8 Tu4 d exc®llent  agreement  with  the  numerical  results  ob- 
tained. The  theory,  however,  only  partially  considers  the  influence 
of  gravity  and  neglects  the  effect  of  drag  and  roll  damping  momen? 

In.thls  report,  we  will  show  that  the  coordinate  system  transfor- 
ls  ^necessary  and  that  a proper  linearization  in  the  fixed- 

svTtem  req^reS  the  soluti°n  of  a much  simpler  fourth-order 

cTh  ?llows  the  very  easy  inclusion  in  the  theory  of  the  full 
effect  of  gravity  as  well  as  the  effects  of  drag  and  roll  damping. 

Much  more  importantly,  the  requirement  of  high  stability  factor  is 
eliminated  so  that  the  very  important  case  of  a finned  missile  with 
little  or  no  spin  can  be  studied.  Finally,  the  effective  technique  of 
quasi -1  mean zat ion  wi  1 1 be  used  to  derive  the  effect  of  the  truly  non- 
linear part  of  the  d>pp  terms  on  the  frequencies  and  the  damping  rates 
of  the  motion. 


1.  K.H.  Lloyd  and  D.P . Brown , Influence  of  Gravity  and  Applied  Side 
Foraee  on  the  Stability  of  a Spinning  Projectile,"  Weapons  Reaearah 
EstabUshment  TR  1906 (W),  South  Australia,  November  1977 

tD  *°.536,48-  .JSee  «l80  "Instability  of  Spinning  Projectiles  During 
lerrmnal  Guidance,  Journal  of  Guidance  and  Control  2,  Jan-Feb  1979. 
pp • 65-70.  ) — - 

*FAXldlla"e  ax.es.^  * pitch  and  yaw  with  the  missile  but  roll  so 
that  the  y-axts  vs  always  in  the  horizontal  plane. 
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II.  EQUATIONS  OF  MOTION 


The  plana  trajectory  of  a particle  flying  at  velocity  V and 
trajectory  angle  6^  with  respect  to  the  horizontal  can  be  described 

by  the  equations2 


where 


V' 

V 

■ - CD  - * 

eT 

" ■ ( t v" 

C* 

-£S£C 

2m  lD 

2* 

■ g l V-2  ! 

(2.1) 

(2.2) 


and  where  derivatives  are  with  respect  to  the  nondimensional  arc- 
length,  s.  These  equations  are  good  approximations  for  the  actual 
variation  of  V and  for  a symmetric  missile.  The  moments  and 

transverse  forces  that  have  a measurable  effect  on  the  missile's 
motion  are  usually  expressed  in  missile-fixed  coordinates3  as 


M 


X 


My  ♦ i Mz 


pstv2h% 

p s . *[(♦<  C% 


K 


+ i 


(2.3) 


(2.4) 


2.  C.H.  Murphy,  "Gravity -Induced  Angular  Motion  of  a Spinning 
Mieeile,"  Ballistic  Research  Laboratories  Report  No.  1546 , July 
1971,  AD  730641,  (See  also  Journal  of  Spaoeoraft  and  Rockets  8 , 
August  1971,  pp,  824-828.) 

3.  C.H.  Murphy , " Free  Flight  Motion  of  Symmetric  Missiles ,"  Ballistic 
Research  Laboratories  Report  No.  1216,  July  1963 , AD  442757. 
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(2.5) 


P¥.iF2.-(i)pSV*  v 
where  <p/  = p fc  V 1 

C = ( v + i w)  v'1 

p = (q  + i r)  i V"1 

The  complex  variable  £ locates  the  plane  of  the  velocity  vector  and 
has  a magnitude  that  is  the  sine  of  the  total  angle  of  attack. 

The  roll  equation  can  be  obtained  for  the  roll  moment  of 
Equation  (2.3)  and  differs  from  the  usual  roll  equation3  by  a gravity 
term  that  acts  on  the  dynamic  pressure: 


where 


P"  = (Dx  ♦ g*)  ♦'  + D: 


D.  = 


CD  + k;2  Cl 


Do  = 


k’2  C* 

3 r ^ i 


k = [I  /m  *2]* 


(2.6) 


In  addition  to  the  aerodynamic  force  and  moment  and  the  gravity 
force,  we  assume  constant  amplitude  control  forces  and  moments  that  are 
perpendicular  to  the  projectile’s  axis  and  either  in  the  horizontal 
plane  or  the  vertical  plane  containing  the  missile's  axis.  These 
control  forces,  which  could  be  produced  by  roll-stabilized  canards, 
make  fixed-plane  coordinates  most  suitable  for  the  analysis.  Since 
fixed-plane  axes  pitch  and  yaw  with  the  projectile  but  roll  so  that 
the  y axis  is  always  horizontal,  it  can  be  shown  that  this  system  has 
an  angular  velocity  vector  Hpp  with  components'4: 

ffFP  - ($FP.  r)  (2.7) 


4.  C.H.  Murphy,  " Measurement  of  Nonlinear  Forces  and  Momenta  by 
Meana  of  Free  Flight  Tests  f’’  Ballistic  Research  Laboratories 
Report  No.  974,  February  1956,  AD  93521. 
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where 


FP 


r tan  9 


0 is  the  angle  between  the  missile's  axis 
and  the  horizontal  (S  = q) ; 

A A 

q,  r are  the  pitch  and  yaw  rates  in  the 
fixed-plane  coordinates 


For  this  coordinate  system,  the  derivatives  of  the  linear  and 
angular  momentum  can  be  computed  in  the  usual  way3*4  and  set  equal  to 
the  sum  of  the  external  forces  and  moments.  The  equations  for  the 
transverse  components  can  then  be  given  in  the  form  of  two  first-order 
complex  differential  equations: 


£'  ■ i V ii  = • y cj  £ + (i  cos  0 + l sin  0^)g  l V-2 

“ (2.8) 

+ (FYC  + 1 FZC>  * V2)_1  - 1 ♦pp  5 

* -ip*  = *;2[V  c -ic*i i 

L p«  aJ 

* fc2  cMq  * c;  * s*]  s <2-9> 

- 1 ki2cM.  * 1 4p  «> 

a 

* k;2  (V  * 1 Mzc>("  v2>'‘  - 1 ♦pp  i 

where 

P * I <t>'/I 
x y 

y ■ u V 1 = the  cosine  of  the  total  angle 
of  attack 
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The  starred  coefficients  are  of  the  order  10  while  the 
dimensionless  control  forces  and  moments  will  be  limited  to  at  most 
lO-4  . Thus,  products  of  these  terms  can  be  neglected  when  p is 
eliminated  between  Equations  (2.8)  and  (2.9)  to  obtain: 


C" 


+ [H  - g* 


- - i 


j 

PH 


(2.10) 


- [M+iPTH=4>  + ci  + C 


where 


2 i 


TP 


M 


KFP 


CP  - *Fp)  + i 


[*FP  (H 


-*) 

Y 


+ <f> 


pp]  } 5 


G = [P  cos  0 - i p sin  0]  g £ V 2 + g*  (£'  - i P £) 


C = i [y  k'2  (Myc  + i Mzc)  - (P  - i y'y'V(Fyc  + i FZC)J  (mV2)’1 
= i Y k’2  (Myc  + i Mzc)(m  V2)'1 


III.  LINEARIZED  SOLUTION  FOR  CONTROL  MOMENT 

For  simplicity,  we  will  first  neglect  the  gravity  terms  in 
Equation  (2.10)  and  consider  the  linear  approximation  to  ♦ . In  doing 
this,  it  is  most  important  to  remember  that  Equation  (2.10)  predicts  a 
steady-state  equilibrium  angle 
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K • s.  * 1 ' htt  ?t  (3-1) 

If  the  small  force  terms  are  neglected,  the  real  part  of  Equation 
(2.8)  yields 


Y r A V 1 ■ - 8'  ♦ *'Fp  a (3.2) 

Equation  (3.2)  is  now  multiplied  by  tan  6 and  solved  for  4>pp  . 

6 - (3.3) 


Now 


Y + a tan  6 


0 ■ / q dt  ■ e ♦ (a  - o.)y-1  (3.4) 

© © 


♦J.p  « a { tan  0e  + y^1  («  * «e)  ♦ Y^1  (a  §e  tan  0e)  (B  - §e) 
♦ aj  (a  - ae)2  + a2  (a  - (0  - 6e)  + a 

sag'  tan  0 

C 

where 


(6  - 6/} 
(3.5) 


a = [y  + a tan  0 1 
1 e e e 


-1 


a,  - [2Se  - (1  - 8j)  a tan  ©e3(2Y|)_1 


a2  * (1  + a Ye)6e  Ye 


-3 


a3  - (1  - 52  + 2a  Ye  8J)(a  tan  0ft)(2  y*)'1 


Since  a linear  analysis^is  concerned  with  small  amplitude  motion 
about  the  equilibrium  angle  5 , ♦ is  expanded  in  powers  of  C - £ 

and  its  derivatives.  The  linear  part  of  this  expansion  is: 


* • i a tan  0fi  [l"  * (H  - i P)8']Ce  (3.6) 


The  usual  solution  to  the  linearized  Equation  (2.10)  neglecting  $ is 


where 


A A 


i<f>5 


l ^2 

' - 5e  * ‘l  * * Ka  • 


<p'.  = (1/2)  (P  ± /P2  - 4M) 


K'./K.  = X.  = 
J J 3 


_J !l 

2 - P 


(3.7) 


* 2 4 

It  is  important  to  note  that  |<)/.|  ~ 10’  and  |x^|  ~ 10" 

For  ascending  or  descending  flight,  $ introduces  terms  in 

§'  = U'  + l')/2  and  l"  (£"  + \")/2  . 


The  linearized  Equation  (2.10)  becomes 


1 + (a  i Cg/2)  tan  9e 


][** 


=']■ 


+ (H  - i P)  ('  • (H  + i P T)C 


C - (a  i 


i le/2)  tan  8e  [$"  + (H  - i P)  f'] 


(3.8) 


As  is  shown  in  Reference  5,  the  effect  of  the  conjugate  terms 
in  Equation  (3.8)  is  to  add  two  additional  modes  in  - and  - ^ . 

For  reasonable  values  of  $ , the  amplitude  of  these  modes  will  be 

much  less  than  Kj  and  K2>  If  we  approximate  the  actual  solution  by 

the  two-mode  solution  of  Equation  (3.7)  and  substitute  in  Equation 
(3.8),  the  conjugate  terms  have  no  contribution  to  the  damping  or 


5.  C.H.  Murphy,  "Angular  Motion  of  Spinning  Almost  Syrmetria 

Missiles ,"  U.S.  Army  Armament  Research  and  Development  Command, 
Ballistic  Research  Laboratory  Technical  Report  02121,  November  1978 
AD  A0635S8. 
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frequency  equations  obtained  from  the  coefficients  of  exp  (i  #j)  • 

For  simplicity,  we  make  the  usual  size  assumptions  f Ujl  <<:  I I 1 
| H } <<  , | T | <<  Ujl^  and  retain  only  terms  linear  in  ae  and  Bg  • 

(<fr')2  - P + M ♦ tan  0 [Mo  + P T 8 ]/2  - 0 (3.9) 

j j e e e 

- H ♦'  + P T - *•  . tan  0 [Mg  ■ P T a ]/2 

A * l l * 2 L_  (3.10) 

2 - P 


IV.  LINEARIZED  SOLUTION  FOR  GRAVITY 

For  no  control  forces  (C  ■ 0) , Equation  (2.10)  predicts  a steady- 
state  equilibrium  angle 


l , 

^e  M + i P(T  - g*) 


G/M 


(4.1) 


where 


A 

Since  a 

e 


G = P g £ V-2  cos  9 . 

6 e 

A 

is  zero,  0_  = 6 and  the  linearized  G and  $ become 
i e 

G = G - i g*  P 8 
*=i  g*  (P/M) [8"  + (H  - i P) 8; ] 


Equation  (2.10)  reduces  co 

V * (H  - g*  - i P)V  - (M  + i P T)( 


(4.2) 

(4.3) 

(4.4) 


* G ♦ i g*  (P/M)  [l"  + (H  - i P)  0'  - MB] 


St-  0^- 


The  solution  to  liquation  (4.4)  can  be  approximate*)  by  the  two-mode 
liquation  (3.7)  with  the  result  that  the  second  term  on  the  right  of 
liquation  (4.4)  has  no  measurable  contribution  to  the  frequency  or 
damping. 


I>j)  - P ^ ♦ M - 0 


(4.5) 


(H  - g*)  + P T 

2 ♦'  - P ~ 


(4.6) 


Pifferentiating  liquation  (4.5), 


r' 


(4.7) 


Since  D in  Equation  (2.6)  is  zero  for  a body  of  revolution,  the 
damping  rate  for  a shell  becomes 


■ H M P T [Dj  P ♦ 2 (P  - ♦')  g*  - M'] 


2 - P 

3 


(2  - P): 


(4.8) 


It  is  interesting  to  note  that  as  the  gyroscopic  stability  increases, 
4/  goes  to  P and  goes  to  0 and  the  contribution  of  gravity  to 


' 1 0 ~ 2 ” 

damping  decays  to  zero 


Equation  (4.6)  is  precisely  the  same  as  that  in  Reference  2. 

The  derivation  given  in  Reference  2,  however,  separately  neglected 
the  second  term  in  Equation  (4.2)  and  all  of  $ . The  correct  deriva- 
tion combines  these  terms  and  shows  that  their  combined  effect  can  be 
neglected . 


V.  QUASILINEAR  ANALYSIS  OF  * 


For  horizontal  flight  without  control  forces  or  gravity,  $ is 
cubic  and  was  shown  in  Reference  4 to  cause  a change  in  frequency. 
Numerical  calculations  by  Clark  and  Hodapp6  showed  that  this  frequency 
shift  was  very  well  predicted  by  the  quasilinear  analysis.  In  this 
section,  we  will  derive  the  quasilinear  prediction  for  cubic  4»  and 
ascending  or  descending  flight  with  control  forces  but  no  gravity. 

The  cubic  part  of  the  nonlinear  term  in  y1  on  the  left  side  of 
Equation  (2.10)  can  be  easily  computed: 


(y'/y)  V 


2 


*-i  a h1  s' 


(5.1) 


The  quasilinear  technique  then  assumes  an  undamped  motion  of  the  form 
of  Equation  (3.7)  and  seeks  the  average  in-phase  and  out-of-phase 
contributions  to  the  coefficient  of  exp(i  4>^)  t>y  the  relation 

i SW  - i ♦. 

[F]i  av  = ITT"  / F e ds 

J.av  Q 

If  there  is  no  equilibrium  angle,  Sw  is  2tt  (<{>'  - iji^)"1  » the  wave- 
length of  £ exp(-iij>j)  . For  nonzero  equilibrium  angles,  the  integrand 
has  several  wavelengths  present  and  is  taken  to  be  large  compared 
to  the  largest  of  these. 


Now 

[i  iv 
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i ♦ 


e 


e 


- i 


(5.3) 


6.  E.L.  Clark,  Jr.,  and  A.E.  Hodapp,  "An  Improved  Teohnique  for 
Determining  Missile  Boll  Bate  with  the  Epiayalio  Theory, " Sandia 
Laboratories  SC-DC-70-4768,  April  1970. 


Therefore, 


tv'  r‘  n,.av  ' (1/2>  ^ - ♦,2)  KJ 


+ ($')  K2  6 K“  e 
2 C 1 


-1  ^ *2 


where 


[ iav  ■ r / n a, 
aV  0 


For  a special  value  of  spin,  2 if/  - <j/  and  the  last  term  in  Equation 

(5.4)  has  a nonzero  contribution.  If  we  assume  that  this  special 
value  of  spin  does  not  occur, 


lY  Y av  * <♦',  - ♦',)  *l'2 


[t'y-1  5'1 


2 a v 


- (*'  - *')  K2/2 


’vnrt^  ■ 

y.  , 

• %\  . 1 


, . '%jK| 


A similar  calculation  can  be  done  for  the  nonlinear  part  of  ♦ . 
The  combined  nonlinear  contributions  to  the  first  mode  have  the  form*: 


[y7  y’1  C7  ♦ * - ♦«] 

* l,av 


b,  + (terms  that  are  zero  when  a » £ ■ 0) 

• © © 

(5.7) 


where 

b. 


£y7  y"1  57J  * a2  tan2  0e  £(87)2  (£  - 

- 1-Li  f t't'H  - i,)] 

® L Jl.av 

- ~r e7  (s  - ae)  (i  - $e)i  - —■  r e7/  (s  - oe)d 

e «■  Jl.av  6 


- t 


■l. 


av 


* Cl/25  [(*')2  K2  * P ♦'  K|  ] a2  tan2  «e 

- d y,)'1  (♦;  - ♦')  [a  k>  . (a  - „ej  ♦;  K2] 


A similar  expression  for  the  second  mode  can  be  obtained  by  inter- 
changing subscripts  1 and  2. 


These  average  nonlinear  contributions  can  now  be  added  to  the  co- 
efficients of  exp(i$j)  in  the  usual  derivation  of  the  linear  damping 

rates  and  frequencies.  For  the  special  case  of  no  control  forces  and 
moments,  the  damping  rates  are  unaffected  and  the  frequency  equations 
become 


P *i  + Ml.av  + jtan2  6e  W2  K?  + P *'t  Kf] 

♦ (*'2  - *j)  *7  “ljj*  0 


(5.8) 


*For  eimplioityt  the  quite  email  nonlinear  terms  in  H have  been 
omitted. 
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Table  I.  Assumed  Parameters  for  105mm  Shell 


l - .105  m 

S - .0087  m2 

m • 15  kg 

I * .023  kg-m2 

x * 

Iy  3 -22  kg-m2 

p ■ 1.05  kg/m3 

V * 250  m/s 

p = 1050  rad/s 


0.13 


B 1.7 


C - - 0.012 

*P 


CM  ‘ 3'8 

a 


0.2 


+ C, 


M. 


a 


8 


and  a similar  equation  for  the  other  mode.  These  frequency  equations, 
for  6 ■ 0 and  an  appropriate  choice  for  M,  were  the  equations  that 

gave  the  excellent  agreement  with  the  numerical  calculations  of 
Reference  6. 


VI . DISCUSSION 


According  to  our  analysis,  the  effect  of  gravity  on  damping  rates 
is  to  replace  H by  H-g*  . Equation  (3-10)  for  the  damping  rates  be- 
comes 


- (H  - g*)  ♦'  f P T - ♦j  - (tan  e^M 


P T 


tAt 


2 


(6.1) 


Since  jp  T|  is  usually  much  smaller  than  | M J , the  effect  of  a on  the 

damping  rates  can  be  neglected.  Equation  (6.1)  can  therefore  b§  used 
to  derive  stability  boundaries  for  the  maximum  trim  angles  for  8 . 

For  a gyroscopically  stable  missile  with  positive  spin 

(♦'  > P/2  > ♦'  > 0), 


where 


B < 8 tan  0 < B 

l e e 2 


(6.2) 


Bj  - - (2/M)  £ (H  - g*)  ♦'  - P T + J 


Table  I gives  the  various  parameters  for  a 105mm  shell  and 
Figure  1 gives  the  boundaries  Bj  and  B2  in  degrees  as  functions  of 

the  gyroscopic  stability  factor,  s^  . The  Lloyd  and  Brown  results 

were  limited  to  large  stability  factors  and  their  numerical  calcula- 
tions were  for  s ■ 5.8  . As  we  can  see  from  Figure  1,  the  allowable 
* 8 

range  for  8 becomes  quite  small  for  stability  factors  less  than  2 . 
c 
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Figure  1.  B-  versus  s for 


A particularly  interesting  case  is  that  of  a statically  stable 
missile  with  little  or  no  spin  (M  < 0,  |s  | « 1)  . For  this  case, 

o 


1 . -rf/M 
M 2 ^ 


(6.3) 


For  an  exponential  air  density  and  constant  static  moment  coefficient > 


M'  /M  ® p'  /p  = - o £ sin  0^ 


(6.4) 


where 


P = PQ  e 

Since  = i , the  stability  bounds  reduce  to  the  very  special 
condition: 

|§e  tan  6e|  < B (6.5) 

where 


2(H  - g*)  - o l sin  0T 

”7^ 


It  is  interesting  to  note  that  the  value  of  8 for  the  Copperhead 
missile  is  estimated  to  be  12  degrees. 
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drag  force 
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roll  damping  moment  coefficient 
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Magnus  moment 
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sum  of  the  damping  moments 
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static  moment 
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Y C.  ♦ Cn  , the  normal  force  coefficient 
a D 

cs  * < 


k;2  sf  ci. 


V Fz 


F F 
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G 

G 


g' 

H 


LIST  OF  SYMBOLS 
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transverse  missile-fixed  components  of  the  aerodynamic 
force 

transverse  fixed-plane  components  of  the  control  force 
P g l V"2  cos  0e 

that  part  of  the  fixed-plane  complex  yaw  forcing 
function  due  to  gravity 

magnitude  of  the  gravity  acceleration 

g l V'2  sin  0T 


C'l  - CD  - k't2  k 

a \ q o/ 


J l 

M 


axial  and  transverse  moments  of  inertia 
magnitude  of  the  j-th  modal  arm,  j - 1,  2 
(Ix/m 

(Iy/m  l2)h 
reference  length 


Ykt  CM 


Mx>  My,  Mz  missile-fixed  components  of  the  aerodynamic  moment 


myc*  mzc 


transverse  fixed-plane  components  of  the  control  moment 


m 

P 


mass 

<vy  *' 


P.  H.  r 


missile  spin,  pitch  and  yaw  rates  measured  in  the 
missile-fixed  system 


•"*> 


LIST  OF  SYMBOLS 
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sw 

s 

S 

g 

T 


t 

U,  V,  w 

V 
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Y 

6 

c 


0 


6 


T 


X. 

J 

v 
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missile  pitch  and  yaw  rates  measured  in  the  fixed-plane 
system 

reference  area 

integration  interval  (calibers)  in  an  averaging  process 


nondimensional  arclength  along  the  trajectory 
gyroscopic  stability  factor 


YC*  +Y 
a 


time 

missile-fixed  components  of  the  velocity 
magnitude  of  the  velocity 

fixed-plane  axes,  the  x-axis  along  the  missile's 
longitudinal  axis  and  the  y-axis  always  in  the 
horizontal  plane 

angles  of  attack  and  sideslip  in  the  fixed-plane  system 
u V'1 


fin  cant  angle 


angle  between  the  missile's  axis  and  the  horizontal 
trajectory  angle 


K'j/Kj  , .1=1,2 

(q  + i r)  £ V-1 
(v  + i w)  V-1 


P 


air  density 
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air  density  at  sea  level 

a 

1/(6700  m) 

4 

that  part  of  the  fixed-plane  complex  yaw  forcing 
function  due  to  the  spin  of  the  fixed-plane  system 
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